Propagation of squeezed radiation through amphfying or absorbing random media 
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We analyse how nonclassical features of squeezed radiation (in particular the sub-Poissonian noise) 
are degraded when it is transmitted through an amplifying or absorbing medium with randomly 
located scattering centra. Both the cases of direct photodetection and of homodyne detection are 
considered. Explicit results are obtained for the dependence of the Fano factor (the ratio of the 
noise power and the mean current) on the degree of squeezing of the incident state, on the length 
and the mean free path of the medium, the temperature, and on the absorption or amplification 
rate. 
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I. INTRODUCTION 

Squeezed radiation is in a state in which one of the 
quadratures of the electric field fluctuates less than the 
other [|]j|]. Such a nonclassical state is useful, because 
the fluctuations in the photon flux can be reduced below 
that of a Poisson process — at the expense of enhanced 
fluctuations in the phase. Sub-Poissonian noise is a deli- 
cate feature of the radiation, it is easily destroyed by the 
interaction with an absorbing or amplifying medium [^ . 
The noise from spontaneous emission events is responsi- 
ble for the degradation of the squeezing. 

Because of the fundamental and practical importance, 
there exists a considerable literature on the propagation 
of squeezed and other nonclassical states of light through 
absorbing or amplifying media. We cite some of the most 
recent papers on this topic M-nffl. The main simplifica- 
tion of these investigations is the restriction to systems 
in which the scattering is one-dimensional, such as paral- 
lel dielectric layers. Each propagating mode can then be 
treated separately from any other mode. It is the pur- 
pose of the present paper to remove this restriction, by 
presenting a general theory for three-dimensional scatter- 
ing, and to apply it to a medium with randomly located 
scattering centra. 

Our work builds on a previous paper |1^ , in which we 
considered the propagation of a coherent state through 
such a random medium. Physically, the problem con- 
sidered here is different because a coherent state has 
Poisson noise, so that the specific nonclassical features 
of squeezed radiation do not arise in Ref. ||ll[. Tech- 
nically, the difference is that a squeezed state, as most 
other nonclassical states, lacks a diagonal representation 
in terms of coherent states [|l|,|[. We cannot therefore 
directly extend the theory of Ref. |l|] to the propaga- 
tion of squeezed states. The basic idea of our approach 
remains the same: The photodetection statistics of the 
transmitted radiation is related to that of the incident ra- 
diation by means of the scattering matrix of the medium. 
The method of random-matrix theory [|l2[ is then used 
to evaluate the noise properties of the transmitted radi- 
ation, averaged over an ensemble of random media with 



different positions of the scatterers. 

The outline of this paper is as follows. In Sec. || we first 
summarise the scattering formalism, and then show how 
the characteristic function of the state of the transmit- 
ted radiation can be obtained from that of the incident 
state. This allows us to compute the photocount statis- 
tics as measured in direct detection (Sec. [II), and in 
homodyne photodetection measurements (Sec. IV). The 
expressions in Sees. pi| pA^ are generally valid for any in- 
cident state. In Sec. ^pwe specialise to the case that 
the incident radiation is in an ideal squeezed state (also 
known as a squeezed state of minimal uncertainty, or as a 
two-photon coherent state ||l|,|[). The statistics of direct 
and homodyne measurements are expressed in terms of 
the degree of squeezing of the incident state. The Fano 
factor, introduced in Sec. VI, quantifies the degree to 



which the squeezing has been destroyed by the propaga- 
tion through an amplifying or absorbing medium. The 
ensemble average of the Fano factor is then computed 
usin g ran dom-matrix theory in Sec. VII. We conclude in 



Sec. VIII 



II. SCATTERING FORMULATION 

We consider an amplifying or absorbing disordered 
medium embedded in a waveguide that supports 7V(w) 
propagating modes at frequency w. The conceptual ad- 
vantage of embedding the medium in a waveguide is that 
we can give a scattering formulation in terms of a finite- 
dimensional matrix. The outgoing radiation in mode n is 
described by an annihilation operator a°"*(a;), using the 
convention that modes 1, 2, . . . , A^ are on the left-hand- 
side of the medium and modes TV -I- 1, ... , 2N are on the 
right-hand-side. The vector a°"' consists of the operators 
a^"' , a™* , . . . , 02^ . Similarly, we define a vector a"^ for 
incoming radiation. 

These two sets of operators each satisfy the bosonic 
commutation relations 

[a„(w), aj„(w')] = (5„,„(5(w - u'), [a„(w), a„(w')] = . 

(2.1) 



They are related by the input-output relations |l^-[l5[| 

a°"*(tj) = 5(w)a'"(w) + Q{iu)b{uj) , (2.2a) 

a°"*(w) = S'(w)a'"(w) + V^(tj)ct(u;) , (2.2b) 

where the first equation is for an absorbing medium and 
the second for an amplifying medium. We have intro- 
duced the 2N X 27V scattering matrix S, the 27V x 2N 
matrices Q and V, and the vectors b and c of 2N bosonic 
operators. The scattering matrix can be decomposed into 
four N X N reflection and transmission matrices, 



5": 



r' t' 
t r 



t^.r 



(2.3) 



and 



Reciprocity imposes the conditions t' 
r' = r'^. 

The operators b and c account for spontaneous emis- 
sion in the medium. They satisfy the bosonic commuta- 
tion relations (pTl), hence 



)t = a - ss\ vv^ = 5'S't - a . 

Their expectation values are 

(6]j(w)6m(w')) = Sn„iS{u - uj')f{uj,T) , 

The Bose-Einstein function 

/(w,T) = [exp(nw/fcr) - 1]^^ 



(2.4) 



(2.5a) 
(2.5b) 



(2.6) 



is evaluated at positive temperature T for an absorbing 
medium and at negative temperature for an amplifying 
medium. 

It is convenient to discretise the frequency in infinites- 
imally small steps of A, so that ujp = pA, and treat the 
frequency index p as a separate vector index (in addition 
to the mode index n). For example, a™' = a™*(wp) and 

^np^n'p' — >Jnn' \^p)^pp' • 

The state of the outgoing radiation is described by the 
characteristic function 

Xout(r?) = (: exp [A'/'{a°"'^^ - 7y*a°-*)] :) , (2.7) 

where (: • • • :) indicates the expectation value of a nor- 
mally ordered product of operators a°"* and 0°"*''' (cre- 
ation operators to the left of the annihilation operators). 
The vector rj has elements rjnp- The density operator of 
the outgoing radiation is uniquely defined by the char- 
acteristic function Xout fl) • Similarly, the incoming state 
has characteristic function 



Xin(r/) = (:exp[Ai/2(a-t,,_^*„in)l .^ 



(2.8) 



The characteristic function of the thermal radiation in- 
side an absorbing medium is given by 



Xabs(r/) = (:exp[Ai/2(6t,,_^*5)j .) 
= exp f - ^ ifnpf{ujp, r)?7„j 

\ n,p 



exp(-?7*/?7) • 
(2.9) 



For an amplifying medium, replacing b by c^ and normal 
ordering by anti-normal ordering, one finds instead 



Xamp(?7) == exp (77*/??) 



(2.10) 



Combination of Eqs. dJ) and (U) with Eqs. ( |2^ 



2.10 ) yields a relationship between the characteristic 
functions of the incoming and outgoing states, 

Xout(r/) = exp (-77*(1 - 55t)/r/) xUS^v) • (2.11) 

This relation holds both for absorbing and amplifying 
medi a, b ecause the difference in sign in the exponent of 
Eqs. (2.9) and (2.10) is cancelled by the difference in sign 
between QQt = 1 - ^-^-t and VV^ = -(a - SS^). 



III. PHOTOCOUNT DISTRIBUTION 

The photocount distribution is the probability P{n) 
that n photons are absorbed by a photodetector within a 
certain time r (see Fig. p. The factorial cumulants kj of 
P{n) (the first two being ki =n and K2 = n{n — l) — n^ 
are most easily obtained from the generating function 

00 ^.- / 00 \ 

no = E ^ = 1^ E(i + ^)"^('^) • (3.1) 

j=l -'' \n=0 / 

The generating function is determined by a normally 
ordered expectation value |lq , [l7[ |, 

»7- 2N 

e^i) = (: e«^ :),]¥= dtY,dnaT'Ht)aT\t) . 

(3.2) 

Here d„ G [0, 1] is the detection efficiency of the n-th 
mode and the time-dependent operators are defined as 

/■oo 

aT{t)^{27r)-'/' cL;e-^-'al^\u;). (3.3) 



Discretising the frequencies as described in Sec. ||, one 
can write 





FIG. 1. Schematic illustration of direct detection: Radia- 
tion is incident on a random medium (shaded). The trans- 
mitted radiation is absorbed by a photodetector. 






^:^^ip-p')taTH^p)aT{^P') 



(3.4) 



This expression can be simplified in the limit t — > oo of 
long counting times, when one can set A = 27r/r and use 



e'^^P-P'^'dt = rSpp, 



(3.5) 



Hence, in the long-time limit the generating function is 
given by 



.Fii) 



(:expUA^d„arn^p)«rK) :) 

\ n,p / 

(: exp(^Aa™*t2?aOut) ;) ^ (3.6) 



where we have defined the matrix of detector efficiencies 

^7ip,n'p' — 0"tiVn.ir' On ri' • 



Comparing Eqs. (2/7) and (3_^) we see that the gener- 
ating function F{S,) can be obtained from the character- 
istic function Xout by convolution with a Gaussian, 



,^(0 _ 



dct(- 



77:^ Mf7Xout(?7)exp f -T?*!? S;j , 



(3.7) 



where J drj is an integration over the real and i maginary 
parts of 7y. We now substitute the relation ( ^.11 ) between 
Xout and Xin, to arrive at a relation between F(^) and Xin: 



.n^) 



1 



det(-C7rX>) 



d-nXin{S'^ri) 



1 



xeK^{-r]*V-^ri-r]*{l-SS'^)f-n\. (3.8) 

The fluctuations in the photocount are partly due en- 
tirely to thermal fluctuations, which would exist even 
without any incident radiation. If we denote by i^th(0 
the generating function of these thermal fluctuations, 
then Eq. (p^) can be written in the form 



no = ^th(e) 



In 



det(7rAf) 



d^l Xin(»7) exp {-rfM ^-q) 



i^th(0 = - lndet[l - iV{t - SS^)f] . 
We have defined the Hermitian matrix 



(3.9) 
(3.10) 



(3.11) 



and we have performed a change of integration variables 
from r] to S'^r] [with Jacobian det(S'5^)]. 

The expression ( 3.10| ) generalises the result of Ref . ||l^] 
to arbitrary detection-efficiency matrix T). Returning to 



a continuous frequency, it can be written as (recall that 

A = 27r/T) 



-Plh(e) = -7^ / d^lndetfl 



-iD[t-S{uj)S\u;)]f(uj,T)) , (3.12) 

where D is a 2N x 2N diagonal matrix containing the 
detection efficiencies d„ on the diagonal (Dnm = dnSnm)- 
The first two factorial cumulants are 



.th 



,th 



duj 

2^ 



f{uj,T)trD[l-S{uj)S^uj)] 



(3.13) 



■ r^f{c.,T)ty{D[l-S{u;)S\u;)])' . (3.14) 



Note that all factorial cumulants depend linearly on the 
detection time r in the long-time limit. 

If only the N modes at one side of the waveguide are 
detected (with equal efficiency d), then dn — ioi 1 < 
n< N and d^ ^ d ioi N + 1 < n < 2N, hence 

Fth (0 = / du In det ( 1 - ^d[l - r (tj)r^ (w) 

27r Jo ^ 

~t{u)t\u)]f{u,T)) , (3.15) 



in agreement with Ref. |l8|| . 

The difference F{^) — Fth(C) contains the noise from 
the incident radiation by itself as well as the excess noise 
due to beating of the incident radiation with the vacuum 
fiuctuations. If the incident radiation is in a coherent 
state, then Xin(^) = exp (a* 77 — 77* a) for some vector a 
(called the displacement vecto r) w ith elements anp — 
On^LUp). Substitution into Eq. (B^) gives the generating 
function 

FiO = ^th(0 - a* Ma , 

^Fti,iO + ir dioa*{Lo)SHLu)x 

27r7o 

X l-^D[l-S{uj)S''{Lu)]f{LU,T)\ DS{iu)a{Lu) 

(3.16) 

The first two factorial cumulants are 
duj 



Kl 







27r 
"^ duj 



a*{uj)S\uj)DS{uj)a{uj) + nf 



th 



(3.17) 



K2 - 2r / ^fiu;,T)a*iu)S^{Lo)D[t - S'(L^)5t(^)] 
Jo 27r 

X DS{uj)a{uj) + K*'' . (3.18) 

If the incident coherent radiation is in a single mo de w p 
and monochromatic with frequency uiq, then Eqs. ( |3.16| - 
3.18) simplify for detection in transmission to 



FIG. 2. Schematic illustration of homodyne detection: At the left, radiation is incident on a random medium (shaded). At 
the right, a strong coherent beam is superimposed onto the transmitted radiation, and the combined radiation is absorbed by 
a photodetector. 



F{0=F,^{0 



^T^dlo 



t^(l-^d[l-rr^ -tt^]f{uJo,T)) t 



Ki = hTd[th],n^^rao + «! 



th 



momo 

(3.19) 
(3.20) 



duj 



nf^rd I --f{uj,T)tY[t~r{uj)r\Lo)~t(uj)t\uj) 



27r 



^2 = 2/ord2/(^o,T) [t\t - rr^ - tt^)t\^ 



(3.21) 

(3.22) 



Kf = Td" 



dhJ 



— f(u;,T)tY[t-r{Lo)r^{u) 



27r 



-t{uj)t\uj)f . (3.23) 

Here /q = (27r)~^ L dw |ap is the incident photon flux 
and the matrices r and t without frequency argument are 
to be evaluated at frequency loq. These are the results of 
Ref. p. 



IV. HOMODYNE DETECTION 



The photocount measurement described in Sec. Ill 
(known as direct detection) cannot distinguish between 
the two quadratures of the electric field. Such phase de- 
pendent information can be retrieved by homodyne de- 
tection, i.e. by superimposing a strong probe beam (de- 
scribed by operators aP™'"') onto the signal beam (see 
Fig. g). The total radiation incident on the detector is 
described by the operator 



^total 



,1/2 out 



(1 



^Nl/2^probc 



(4.1) 



where the factor ^/n accounts for the attenuation of the 
signal beam by the beam splitter that superimposes it 
onto the probe beam. (For simplicity we assume a real 
scalar k, more generally k would be a complex coupling 
matrix.) 

The characteristic function of 0*°*^' is the product of 
the characteristic functions of a°"' and aP'^°'°'^ . We as- 
sume that the probe beam is in t he co herent state with 
displacement vector (3. From Eq. ( |2.11 ) one gets 



Xtotai(?y) =exp 



-Kr/*(l-55^)/77 



The generating function Fhomo{Cj of the photocount dis- 
tribution in homodyne detection is given by [cf. Eq. (^)] 



exp[Fhomo(C)] 



(: exp(^Aa*°*'^it2?atotai) 
exp[^(l-K)/?*I?/?] 



expfAi/V'«(l-'«)^ 



(4.3) 



(4.4) 



In the second approximate equality we have linearised 
the exponent with respect to a°"*, which is justified if 
the probe beam is much stronger than the signal beam. 
The remaining expectation value has the form of a char- 
acteristic function if we take ^ purely imaginary, so that 
i* = -^. The result is 



+ k(1 - k)^^/3*P(1 - SS^)fVl3 

+ lnXin(v/«(l-K)e^^2?/3). (4.5) 

I n the second equation we have substituted the relation 
( ^.llD between Xout and Xin- 



V. SQUEEZED RADIATION 

We consider the case that the incident radiation is in 
the ideal squeezed state \e,a) = CS\0) ||l|,||, obtained 
from the vacuum state |0) by subsequent action of the 
squeezing operator 



S = exp [iA (a'"e*a'" - 
and the displacement operator 

C = exp fAl/2 (flinty 



lint^^inf^j 



a ' ea 



(5.1) 



(5.2) 



As in the previous sections, we have discretised the 
frequency, tUp = pA, and used the vector of opera- 
tors a™p = a^^{ujp). The complex squeezing parameters 

e„(ijj) = /9„(a;)e*'^"*^"-' are contained in the diagonal ma- 
trix e with elements e„p,„'p' = £n{'^p)Snn'Spp'- Similarly, 
the vector a with elements a„p = anitOp) contains the 



+ (1 - k)1/2(/3*^ _ rj*f3) Xin(5'^Ki/\) (4 2) displacement parameters. 



The characteristic function of the incident radiation is 
given by HI] 

Xin(?7) = exp [a*T] - 77*0; - jT]{e~'''^ sink 2 p)r] 

-Iif{e"l'smh2p)r]*-T]*{smh'^p)ri] . (5.3) 



According to Eq. ( 2.11 ), we thus find for the character- 
istic function of the outgoing radiation 

Xontiv) =exp{a*S^'n-7j*Sa- lj]S*{e-"l'smh2p)S^r] 
-lr]*S{e''^ smh2p)S^T]* -v*[f- S{f - sinh^ p)S^]jj) . 

(5.4) 

The generating function F{^) of the photocount distri- 
bution is obtained from Xin by convolution with a Gaus- 
sian, cf. Eq. (13). We find 



i^(C)=i^th(0-5lndetX-i 



a 
a 



X' 



Ma 
M*a* 



(5.5) 



where the matrix X is defined in terms of the matrix M 
by 



X = 1 



M sinh p 
-M*e-"l' cosh p 



(5.6) 



-Me^-^coshp 
M*sinhp 

sinh p 
sinh p 



If squeezing is absent, p — 0, hence X = t and Eq. ( |5.5| ) 
reduces to the resuh ( 3.16 ) for coherent radiation. For a 
squeezed vacuum (a — 0) one has simply F{£_) — Fth{£,)~ 
ilndetX. 

If the radiation is incident only in mode mo, then 
we may co mput e the matrix inverse and the determi- 
nant inEq. (5.5) explicitly. The matrix M{uj) defined in 
Eq. (3.11) may be replaced by its mQ,mQ element. 



M„„,„„(c^) = m = -e S^{t - ^D[t - SS^f) 'DS 



mo mo 

(5.7) 



Note that m is real, since it is the diagonal element of a 
Hermitian matrix. The resulting generating function is 

no = ^th(o 

1 / ClUJ , - 1 2 2 ■ 1 2 \ 

■- ^T I — ln(l + 2?7tsmh p~rn smh p) 
Jo 27r 

dio 2 1 

——m\a\ ; — 2 :r~. — 2 — ^ 

2n 1 -|- 2m sinh p — m^ sinh p 

X I 1 -|- TO sinh p[sinh p + cosh p cos(2 arg a — (/))]) 

(5.8) 

The first two factorial cumulants, for detection in trans- 
mission, are 



/•OC 7 

Ki = K^ + Td — (|ap + sinh^ p)[t^t], 

Jo 27r 



(5.9) 



K2 



.th 



+ 2Tcf 



dhJ 

2^ 



(|a|2 + sinhV)/ 
X [t\l - rr^ - 



tt'')t]r 



rd' 



cLu 

2^' 



[i^ilmomo [|a coshp - a* c'^ sinhpl^ 



|ap + sinh^p(cosh^p + sinh^/9)] , (5.10) 



where nf^ and k^ are given by Eqs. ( |3.21 ) and ( 3.23 ). 

The generati ng f unction for homodyne detection fol- 
lows from Eqs. (|4.5| ) and (p^), 

i^homo(0-e(l-«)r2?/3 

+ ^y/K{l~K){a*S''V/3 + P*VSa) 

- iC^K(l - K)[pVS*{e-''''smh2p)S''Vf3 

+ (3*VS{e''^ sinh 2p)S'^Vp*] 
+ C^k(1 - K)(3*V[f - S{f - sinh^ p)S'']V(3 . (5.11) 

All factorial cumulants except for the first two vanish in 
the strong-probe approximation. We may simplify the 
generating function by assuming that the signal beam is 
incident in a single mode TOq and that the probe beam is 
also in a single mode uq. For detection in transmission 
one then has the factorial cumulants 



Ki — rd 



dio 
2 



vr V 



1 



K2 



-tk(1 — K)d' 
f 2rK(l - K)d'^ 



+ 2 V«(l - K)Re[a(3*tn,rno] ) , (5.12) 
doj , 



2tt 



Re[f3* 



t^ 1 

■'nomoJ 



sinh 2p 



"^1 
2n 



\t„grnof sinh^ p 



rr' 



tt^)r 



(5.13) 



VI. FANO FACTOR 

For the application of these general formulas we focus 
our attention on the Fano factor JF, defined as the ratio 
of the noise power P — r~^varn and the mean current 



T^P/I^l + K2/K1 



(6.1) 



(We have assumed the limit r ^ 00.) For coherent radia- 
tion JF = 1, corresponding to Poisson statistics. Thermal 
radiation has T > 1 (super-Poissonian). Nonclassical 
states, such as squeezed states, can have J- < 1. 

We assume that the radiation is incident in a single 
mode Too and is detected in transmission (equal efficiency 
d per transmitted mode). We consider a frequency- 
resolved measurement, covering a narrow frequency in- 
terval around the central frequency wq of the incident 



radiation. The thermal contributions k'^ and k^ may 
then be neglected, since they are spread out over a wide 
frequency range. The incident radiation has Fano fac- 
tor jFjn , measured in direct detection with unit efficiency. 
For squeezed radiation, one has 



^in = 1 



1 



al 



sinh p 



I a cosh p — a* e^'^ sinh p\^ 



|ap +sinh^p(cosh^p + sinh^p)] . (6.2) 



We seek the Fano factor of the transmitted radiation, 
both for direct detection (.Fdiroct) a nd fo r hom o dyn e de- 
tection (.?"homo)- Combining Eqs. ( |3.9[) and ( 5.1C ), we 
find for direct detection 



'^direct — -L V (J'\T' t Jniomo \*^in ^) 



2d/(c^o,T)i 



{tH)r 



(6.3) 



The first term is due entirely to the incident radiation. It 
is absent for coherent radiation (because then J^i^ — 0). 
The second term is due to the beating of the incident 
radiation with the vacuum fluctuations. It is independent 
of the incident radiation and was studied in detail in 
Rcf. 0. 

The Fano factor in the strong-probe appr oximation 
(|/3| -^ cx)) follows from Eqs. ( ^.12D and ( |5.13| ), with the 
result 

-^homo - 1 = 2dK\tnomo P sinh^ p 

+ 2dnf{iOo,T){t ~ rr^ ~ tt^)n„no 

- dKRe[e^(*-2-g/')4„J sinh2p . (6.4) 

In the strong-probe approx ima tion, it is independent of 
a and |/3|. Similarly to Eq. (6.3), the first term is entirely 
due to the incident radiation, vanishing for coherent ra- 
diation [p = 0), and the second term is due the beating 
with vacuum fluctuations. The additional third term de- 
scribes the effect of the phase of the probe beam on the 
measurement. Typically, in a measurement one would 
vary the phase of the probe beam until the Fano factor 



is minimised, which occurs when arg/3 = -kef) 
The resulting Fano factor ^^^^ is given by 



-argt„ 



.77 nun 
*^honio 



= 1 - 2dK\tno,no\'^e ''sinhp 

+ 2dKf{ujo, r)(]I - rr^ - tt^)nono ■ (6.5) 



VII. ENSEMBLE AVERAGES 



random-matrix theory [|2|. Ensemble averages of mo- 
ments of rr^ and tt^ for iV ^ 1 have been computed by 
Brouwer pC| , as a function of the mean free path I and 
the amplification (absorption) length £^a = y/Dra, where 
I/tq is the amplification (absorption) rate and D ~ cl/i 
is the diffusion constant. It is assumed that both ^a and 
L are small compared to the localisation length Nl but 
large compared to the mean free path I. Obviously, this 
requires a large number N of propagating modes. The 
relative size of L and ^a is arbitrary. 

As sample-to-sam ple fl uctuations are small for A^ ^ 1, 
we can take in Eq. ( |6.3| ) the averages of numerator and 
denominator separately. The dependence on the in- 
dex 777-0 of the incident mode drops out on averaging, 
(• ■ •)momo = N~^{iT: ■ ■ ■). For an absorbing disordered 
waveguide, we find 



'^direct -L 



Aid 



3^a sinh s 
2s + cotanhs 
sinh s 



(^in - 1) 



^/K,r) 

i cotanh s — 1 
sinh s 



sinh s 



(7.1) 



We have abbreviated s = L/^a- In the limit of strong 
absorption, s — *■ cxd, the Fano factor approaches the uni- 
versal limit [0 J^dircct = 1 + irf/- The Fano factor T-m 



is gi ven by Eq. (6.2) for an incident squeezed state, but 
Eq. (7.1) is more generally valid for any state of the in- 
cident radiation. 

The result for an amplifying disordered waveguide fol- 
lows by the replacement Tq — > — r^, hence ^a — *■ iS,a'- 



T, 



direct 



= 1 



Md 



iia sm s 
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(7.2) 



The Fano factor diverges at the laser threshold s = tt. 
The function f{ujo, T) now has to be evaluated at a neg- 
ative temperature. For a complete population inversion 
of the atomic states / ^ —1. 

The minimal Fano factor in homodyne detection is 
given by Eq. (|^). The average (|t„omoP) is again in- 
dependent of the mode indices, hence it can be replaced 
by N~'^{trtt''). For an absorbing waveguide we find 



•'"homo ^ 



8ldK 



+ 



3N£_a sinh s 
8ldK 



e '' sinh p 



5ia 



'■f{^o,T) 



cotanh s 



1 



sinhs 



(7.3) 



The expressions for the Fano factor given in the previ- 
ous section contain the reflection and transmission ma- 
trices of the waveguide. These are iV-dimensional ma- 
trices that depend on the positions of the scatterers in- 
side the waveguide. The distribution of these matrices 
in an ensemble of disordered waveguides is described by 



and for an amplifying waveguide 
8ldK 
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(7.4) 
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FIG. 3. Average Fano facto r .7-dircct for direct detection as a function of the length of the waveguide. The left panel is for an 
amplifying medium [Eq. (7.2), / = — 1], the right panel for an absorbing medium [Eq. ([7.1]), / = 10~^]. In both cases we took 
l/S^a = 0.1, d = 1, and values of jTin increasin g fr om t o 3 in steps of 0.5. The dotted parts of the curves are extrapolations in 



the range L < I that is not covered by Eqs. (7.1) and (7.2) 



iTomo requires 



Measurement of the ensemble average J-^ 
that for every sample the phase of the probe beam is 
re-adjusted so as to minimise the Fano factor. If the 
phase of the probe beam is fixed, the random pha se o f 



U. 



will average to zero the third term in Eq. (3.4) 



' sinhp. 



In Eqs. (7.3) and (7.4) this amounts to the substitution 



A graphical presentation of the results (7.1)-(7.4) is 
given in Figs. and 0. For the absorbing case we have 
taken / = 10~ (corresponding to optical frequencies at 
T = 3000 K). For the amplifying case we have taken 
/ = — 1 (complete population inversion). The formulas 
above cannot be used for L <l. The values of JFdircct and 



•^homo for i = can be read off from Eqs. (BJ) and (3.5) 



^d 



ircct 
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An extrapolation to i = is shown dashed in Fig. ||. 

The common feature of the Fano factors plotted in 
Figs. H and ^ is a convergence as the length of the 
waveguide becomes longer and longer. For an absorbing 
medium the L — > oo limit is independent of the state of 
the incident radiation. For an amplifying medium com- 
plete convergence is pre-empted by the laser threshold at 



VIII. CONCLUSIONS 

In conclusion, we have derived general expressions for 
the photodetection statistics in terms of the scatter- 
ing matrix of the medium through which the radiation 
has propagated. These expressions are particularly well 
suited for evaluation by means of random-matrix theory, 
as we have shown by an explicit example, the propaga- 
tion of squeezed radiation through an amplifying or ab- 
sorbing waveguide. The sub-Poissonian noise that can 
occur in a squeezed state (characterised by a Fano fac- 
tor smaller than unity) is destroyed by thermal fluctua- 



tions in an absorbing medium or by spontaneous emission 
in an amplifying medium. The theory presented here 
describes this interaction of nonclassical radiation with 
matter in a quantitative way, without the restriction to 
one-dimensional scattering of earlier investigations. 
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